:arXiv:1111.3295vl [astro-ph.GA] 14 Nov 2011 



B 



hi 



•a 



I 

a" 
o 



(T) 



Astronomy & Astrophysics manuscript no. ISOGAL-CMZ 


©ESQ 2011 


November 15, 2011 





Recent star formation in tlie inner Galactic Bulge seen by ISOGAL 

II - The Central Molecular Zone 

K. Immer^'^, F. Schuller\ A. Omont^, and K. M. Menten' 

' Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany 

~ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, 02138 Cambridge, MA, USA 

^ Institut d' Astrophysique de Paris, CNRS, 98 bis, Bd Arago, F-75014 Paris, France 

Received 09/08/2011; Accepted 08/1 1/201 1 

Abstract 

We present 5-38 yum spectroscopic observations of a sample of 68 ISOGAL sources with unknown natures, taken with the Spitzer 
Infrared Spectrograph. Based on the characteristics and the slope of their spectra we classified the sources as young or late-type 
evolved objects. These sources were selected to test selection criteria based on the ISOGAL [7]-[15] color and the spatial extent 
parameter tri^. We revised these criteria until they reliably distinguished between young and late-type evolved objects and then 
applied them to all ISOGAL sources in the central molecular zone (CMZ), resulting in the selection of 485 sources believed to be 
young. Furthermore, we added 656 Midcourse Space Experiment (MSX) sources to the CMZ sample that fulfilled Fe/Fo > 2 with 
FD and Fe being the flux densities in the D (15 jum) and E (21 fim) bands. After obtaining jjj conversion factors, we calculated the 
bolometric luminosity, Lboi, values for the CMZ sample and subsequently the masses of the sources. Applying a Kroupa initial mass 
function, we derived the total mass in young objects that has been formed over the last 1 Myr, resulting in an average star formation 
rate of 0.08 solar masses per year for the CMZ. 

Key words. Galaxy: center - Galaxy: stellar content - Stars: formation - H II regions - Infrared: stars 



1. Introduction 

The central few hundred parsecs of the MiUcy Way present 
an exceptional environment in our Galaxy. This region pro- 
duces around 5-10% of the infrared luminosity of the Galaxy 
and contains ca. 10% of its neutral gas (~5-1 0-10^ M0) (e.g. 
iGustenll 1 9891: iRodriguez -Fernandez et al. ||200i. The conditions 
in this central molecular zone (CMZ, the inner + 1.5° x ± 0.5° 
around the Galactic Center) are very extreme, the temperature, 
the velocity dispersion, the pressure and the estimated mag- 
netic field strengths are m uch higher than in the Galactic disk 
(iMorris & Serabvnlll996h . From a statistical point of view, a 
study of the star formation activity in the CMZ can uncover if 
star formation in this region is enhanced or suppressed by these 
extreme conditions. 

The CMZ harbors many sites of massive star formation. The 
prominent emission of giant molecular clouds from infrared to 
radio wavelengths indicates recent bursts of star formation near 
the Galactic Center However, most star formation studies in the 
CMZ are focused on the most massive and highly active star 
clusters near Sgr A* or the giant H II complexes Sgr B2 and Sgr 
C. 

Thus, in order to obtain an extensive overview of the star 
formation activity in the CMZ, a large number of young stellar 
objects has to be identified in this region. A method which is 
applicable on a large number of sources to identify YSO candi- 
dates is the usage of color selection criteria. However, in order 
to be reliable, photometric color criteria require always a spec- 
troscopic confirmation of the nature of objects in a well selected 
sample. 

At a particular time in its evolution a young stellar object 
(YSO) begins to heat its dust in the surrounding birth cloud 



which then re-emits this absorbed energy at infrared wave- 
lengths. Studying massive star formation at infrared wavelengths 
(preferably with spectroscopic methods) gives therefore insight 
into the characteristics of the circumstellar dust as well as the 
properties of the embedded star. 

Large surveys combining several bands in the near- and mid- 
infrared are important tools to carry out systematic studies of 
young stellar populatio ns. In this pu blication, we will use data 
of the ISOGAL (Omont et al.ll2003l) and the Midcourse Space 
Experiment (MSX, Price et al. 2001) surveys. The ISOGAL sur- 
vey comprises observations of the inner Galaxy at 7 and 15 /urn, 
taken with the infrared camera ISOCAM on board the Infrared 
Space Observatory (ISO) with a spatial resolution better than 6". 
The MSX project is a multi-wavelengths infrared survey of the 
entire Galactic plane in four spectral bands at 8, 12, 15, and 21 
jum with a spatial resolution of ~18". 

ISchulIer' et al.l (l2006i) (hereafter: Paper I) defined selection 
criteria based on the [7]-[15] ISOGAL color (where [7] and 
[15] are the magnitudes at 7 and 15 yum, respectively, measured 
with ISO), as well as the magnitude uncertainty at 15 fim, cris 
(residual between the fi tted point spread fu nction and the ac- 
tual source profile, see ISchuIler et alJ l2003h . which serves as 
a parameter for the spatial extension of the source (explained 
below). These criteria were established from a group o f point 
sourc e s from the ISOGAL point source catalogue ( SchuIIe r st al.l 
I2OO3I : lOmont et al. 2003) with known natures, dividing them 
into young objectfland late-type evolved objects (LEOs). Since 



'In this publication, the term 'young object' refers to YSOs (a phase 
in the star formation process before the ionization of the surrounding 
birth cloud) and H II regions (a later phase where the star starts to ionize 
its dust cocoon). 
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young objects are embedded in a dust cocoon, they appear more 
extended at 15 fim than LEOs. As explained in Paper I, a high 
value of cri5 indicates that the source is marginally resolved at 
15 jum. To be present in the PSC, such sources still have angular 
diameters < 30" but values of 10"-20" (FWHM) are typical. 

In this publication, these selection criteria are tested on a 
sample of sources with unknown evolutionary stages. We will 
determine the nature of these sources from the characteristics 
and the slope of their mid-infrared spect ra which were obta ined 
with the Spitzer Infrared Spectrograph jHouck et al.ll2004) . By 
means of this sample of unknown sources we will refine the se- 
lection criteria until they reliably distinguish between young and 
late-type evolved objects. 

The final goal is to apply a high quality version of the ISO 
selection criteria on all ISOGAL sources in the CMZ in order to 
find young objects and to characterize the star formation activity 
near the center of our Galaxy. 

This paper is divided into two parts. The first part comprises 
the analysis of the infrared spectra. In Sections 12.11 12.21 and 
I2.3l we present the selection of sources for our test sample, the 
Spitzer observations of the infrared spectra and the data process- 
ing of the spectra. In the following subsection, the classification 
of the studied sources is described. In the last subsection of this 
part, we will give refined selection criteria, which will then be 
applied to all ISOGAL sources in the CMZ in the second part of 
this publication in order to find young objects within the inner 3° 
of our Galaxy. In Section [372l we derive the bolometric luminosi- 
ties as well as the masses of all sources in our CMZ sample and 
finally determine the average star formation rate in the CMZ. 



2. Verification of tiie ISO color selection criteria 

2. 1 . Source selection 

Infrared colors are a powerful tool to distinguish between 
various classes of objects. Based on fluxes at 12, 25, and 
60 //m measured with the Infrared Astronomy satellite (IRAS), 
[Wood & Churchwell (1989) defined criteria to identify ultra- 
compact HII (UCHII) regions based on far-IR colors. In the 
1990s, mid-IR surveys conducted with MSX and ISO covered 
the Galactic plane with much high er spatial resolution and sen- 
sitivity than IRAS. iLumsden et al.l (120021) defined color criteria 
based on MSX and the near-IR 2MASS survey to select candi- 
date massive YSOs. Similarly, from identifications of ISOGAL 
point sources with known objects, it has been found that mas- 
sive YSOs a nd UCH II regi ons show [7]-[15] colors greater 
than 2 mag (Felli et al.1 120021) . Moreover, most sources with 
[15] < 4 mag and [7]-[15] > 2 mag have MSX counterparts with 
flux ratios F21 ;jm/Fi5 fjm above 2, making their interpretation as 
YSO very likely, since YSOs are characterized by strong rising 
continuum emission in the IR range. In addition, in Paper I it 
has been shown that the spatial extent at 15 fim, as it appears in 
the ISOGAL catalogue, can help in determining whether sources 
brighter than 300 mJy at 15 fj.m (or [15] < 4.5 mag) are young 
objects or late-type evolved stars. Infrared sources with ISOGAL 
characteristics similar to those of known massive young objects, 
but without any known counterpart, are thus good massive YSO 
candidates. 

Not all of the sources in the ISOGAL catalogue have been 
detected at both of the covered wavelengths (7 and 15 yum). Since 
sources brighter than 8 mag at 7 yum should have been detected 
by the ISOGAL survey, a lower limit for the [7] magnitude of 
8 mag was assigned to sources which were not detected at 7 fim. 
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Figure 1 : Footprint of the IRS slits shown on a 24 //m image 
around J 174505.6-29 101 8. In that case, the map consists of 6 
pointings with the SL module, and 2 pointings with the LL mod- 
ule. Positions in the slits where background spectra were ob- 
tained are shown in blue. 



Near the Galactic Center. loiass et al ] (1200 Ih published a cat- 
alogue of long period variable stars, which also appear as bright 
ISOGAL sources. In the 0. 1 deg^ regi on of the ISOGAL survey 
overlapped by the iGlass et aP (1200 ih observations, 95 sources 
brighter than 5.25 mag at 15 fim show [7]-[15] colors larger 
than 1.8 mag. Sixteen of these s ources can be matched with 
long period variable stars from the lGlass et al ](2001)catalo gue, 
and were removed from our sample. After rejecting another six 
sources associated with known OH/IR stars and five showin g 
near-IR spectra typical of AGB stars (ISchultheis et al.l l2003h . 
we finally selected 68 ISOGAL sources as targets for follow- 
up spectroscopic observations with Spitzer. This test sample is 
complete down to a 15 yum flux density of 140 mJy over a 0.1 sq. 
deg. area. 

Despite our effort to remove known LEOs from our test 
sample, we realized later that two of our targets (J174602.8- 
29 0359, J174623.0-285845 ) have been identified as OH/IR stars 
bv'Siouwerm an et"an (Il998h . These sources have been helpful in 
establishing the identification criteria for LEOs in our test sam- 
ple. 

2.2. Observation and Inspection of the Data 

The spectroscopic data were obtained in March 2005 with 
the t wo low resolution modules of the Infrared Spectrograph 
(IRS; 'Houck et al."'2004^), onboard the Spitzer Space Telescope 
(Werner et al. 2004). The 68 sources were observed in the 
IRS-mapping mode, with the low-resolution modules SL 
(Short-wavelength, Low-resolution) and LL (Long-wavelength, 
Low-resolution). Both modules are subdivided in three sub- 
modules that cover different parts of the two wavelength ranges 
(SL2: 5-7.5 fim, SL3: 7.3-8.7 //m, SLl: 7.5-15 jum, LL2: 
14-21 fim, LL3: 19.4-21.7 yum, LLl: 21^0 yum). Each map 
consisted of between 3 and 17 pointings with the SL module, 
and between 1 and 4 pointing(s) with the LL module, depending 
on the size of the source as seen in the ISOGAL image. As 
an example. Fig. [T] shows a footprint of the IRS slits on a 
Spitzer/MIPS 24 yum image around the source J174505.6- 
291018. The ramp duration was set to 6 s for each module. 
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and two cycles were required per ramp. In most cases, this was 
sufficient to achieve S/N ratios of 100 or more over the entire 
wavelength range. 

We began with a first inspection of the data. Depending on 
the detection of the ISOGAL sources at short (5-15 fim) or long 
(15^0 /im) mid-infrared wavelengths in the spectroscopic data, 
we divided the sources into different groups: 

- A) ISOGAL sources, detected as one source in the SL as well 
as in the LL data at the same position (49 sources) 

- B) ISOGAL sources, only detected in the LL data (8 sources) 

- C) ISOGAL sources, detected neither in the SL nor in the LL 
data (4 sources) 

- D) Confusion of several Spitzer sources in the SL or LL data 
at the position of the ISOGAL source (7 sources) 

Only the sources in groups A and B were processed further, 
giving a total of 57 sources. Since the group C ISOGAL sources 
are not detected with Spitzer at higher resolution, these sources 
are probably false detections in the ISOGAL catalogue. 

2.3. Data reduction 

The data reduction was conducted with the Spitzer software 
IRSCLEAN0 (to clean rogue pixels in the data) and SPIC^I (to 
extract the spectra). The spectra are calibrated by SPICE with 
calibration files which are automatically selected by the software 
(for more information about the calibration see the SPICE User's 
Guid^. 

A direct background subtraction cannot be conducted with 
SPICE. Therefore, the source and the background spectra were 
determined separately and the background subtraction was con- 
ducted later. Since the infrared background in the Galactic 
Center region varies greatly on short spatial scales, the back- 
ground spectra were obtained from the same subslit as the source 
spectra. Several background spectra were extracted from differ- 
ent positions and pointings around the source (the number of ob- 
tained background spectra varies from source to source, depend- 
ing on the infrared source density around the target sources). The 
blue lines in Fig. [T] show the positions where background spec- 
tra were obtained around J 174505.6-29 101 8. In order to obtain 
a mean background spectrum at the position of the source, the 
weighted average of the background spectra and the associated 
errors were determined. 

Per subslit and submodule (SL2, SL3, SLl, LL2, LL3, LLl), 
the source spectra were extracted from those pointings where 
the emission of the source could be distinguished from the back- 
ground emission and the emission of other infrared sources. Per 
pointing, the source spectra of the two cycles were averaged and 
the mean background spectrum was subtracted. Per submodule, 
all background subtracted source spectra were added together 
Then, the background subtracted source spectra of the six sub- 
modules were assembled to produce the final source spectrum. 

Seventeen sources of the test sample show flux jumps be- 
tween the SL and the LL module spectra around 14 fim. These 
are due to extended emission that falls in the slits in combination 
with the different slit widths of the two modules covering areas 
of different sizes around the sources. 

Twentyseven of the 57 test sample sources were only de- 
tected at 15 jum in ISOGAL and a lower limit of 8 mag was as- 



signed to their [7] magnitudes. Therefore, their [7]-[15] colors 
are also only lower limits. Furthermore, it should be mentioned 
that the ISOGAL magnitudes at 7 and 15 //m were not observed 
at the same time for 63% of the 57 sources (sources with right 
ascension, RA, < 17:45:40). Due to possible variability of the 
sources, this could lead to errors in the [7]-[15] parameter 

2.4. Spectroscopic Properties of tlie Test Sample 

The most common features in the IRS wavelength range, de- 
tected in our sources, are silicate absorption bands, emission fea- 
tures from polycyclic aromatic hydrocarbons, and fine structure 
lines from ionized species. A broad silicate absorption band at 
around 9.8 fim is present in all the spectra but the central wave- 
length changes marginally from source to source. Further sili- 
cate absorption features are detected at 18, 22.3 and/or 23 fim 
(e.g. iMolster et alJl2002h . These silicate features are indicative 
of O-rich envelopes with large mass infall rates. 

In a few sources, absorption bands at around 5.95, 6.7, an d 
7.4 yum seem to be present. According to Gerakines e t al.l (Il995h . 
the absorption at 6 yum is partly due to a bending mode of water 
ice. The 6.7 jim absorption ba nd can partly be exp lained by ab- 
sorption through methanol ice (lOuirico et alJl2000l) whereas the 
band at 7.4 /im is possi bly caused by acetaldehyde (CH3HCO) 
ice (Sch utteetalll 998'). 

Emission from Polycyclic Aromatic Hydrocarbons (PAHs) 
requires a strong ultraviolet field to be excited and thus is fre- 
quently found at the edges of H II regions. Eleven of the sources 
show unambiguous emission of PAHs at one or more of the 
following wavelengths: 6.2 fj.m, 7.7 fj.m, 8.6 fim, 11.3 jum, and 
16.4 /urn. However, due to confusion of the ice absorption bands 
at wavelengths shorter than 8 fim, in several sources the PAH 
emission features in this wavelength range could not be unam- 
biguously identified. 

In our source spectra, the detection of PAH emission bands 
is always seen in combination with at least one of the fol- 
lowing forbidden fine structure lines: [Aril] (6.97 jum), [Nell] 
(12.83 yum) and [SIII] (18.7 yum). The most frequent emission 
line, [Nell], is detected in the direction of 21 sources. In five 
of the sources the detection is weak and could be spurious. The 
[Aril] and the [S III] lines are observed in the direction of seven 
and fourteen sources, respectively. Although the occurrence of 
these emission lines often indicates the existence of H II regions 
and therefore the presence of hot massive stars, this emission 
can also arise from the warm, ionized, interstellar medium sur- 
rounding the sources. This fact complicates the background sub- 
traction. The detection of the emission lines should therefore be 
considered with caution. Nevertheless, the lines can still be used 
as a diagnostic tool. 

In several sources, the S(l) and S(3) rotational emission lines 
of H2 at 17.03 fim and 9.66 jum (within the broad silicate absorp- 
tion feat ure) are detected which are tracers for w arm molecular 
gas (e.g. lKaufman et alJl2006trFiolet et alJlIoTol) . 



: //ss c ■ spitzer ■caltech.edu/ar chanaly/contributed/irs clean/IRSCLEAM_MASK.html| 
■^http ://ssc ■ spitzer ■ caltech.edu/ dataanalysi stools/tool s/sp ice/| 
" ^ihttp : // ssc . spitzer . caltech.edu/postbcd/doc/spice . pdf. 
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Table 1: Detected spectral features of the sources. The first column contains the source names. Columns 2 through 8 give the 
different detected spectral features: PAH emission, emission of the forbidden fine structure lines [Aril], [NellJ, and [SIIIJ, emission 
of the rotational lines S(l) and S(3) of H2, and silicate absorption at 22.3 and/or 23 jum. The last column lists the detection of ice 
absorption bands in the spectra. An unclear identification of a spectral feature was indicated with a question mark. 
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Figure 2: Background subtracted spectrum of the source 
J174505.6-291018. The data of the different subslits are pre- 
sented by the three different colors (green: SL2 or LL2, blue: 
SL3 or LL3, red: SLl or LLl). Denoted are the silicate absorp- 
tion features, the PAH emission bands and the forbidden fine 
structure lines [Aril], [Ne II], and [S III]. 



Table [T] shows the spectral features that were identified in 
the source spectra. The first column contains the source names. 
Columns 2-5 list the detection of PAH emission features as well 
as forbidden fine structure lines ([Aril], [Nell], [SIII]). Where 
the identification of a feature was not clear, question marks were 
set in the table. Columns 6-8 mark the detection of emission of 
the S(l) and S(3) rotational lines of H2 at 17.03 jum and 9.66 fim 
as well as silicate absorption at 22.3 or 23 fim. In the last col- 
umn, the detection of ice absorption bands at around 5.95, 6.7, 
or 7.4 jjm are mentioned. 

As an example. Figure |2] shows the background-subtracted 
spectrum of the source J174505. 6-291018. The three different 
colors represent the data of the different subslits (green: SL2 or 
LL2, blue: SL3 or LL3, red: SLl or LLl). The spectrum shows 
emission lines of [Aril], [Nell], and [SIII] as well as two sili- 
cate absorption features at 9.8 and 18 //m. Furthermore, the de- 
tected PAH emission features are indicated. The spectra of the 
remaining 56 sources are shown in Fig. [3] (online material). 

2.4.1 . Classification of tine 57 sources 

The vast majority of ISOGAL sources brighter than 5 mag at 
15 /zm are either of young stellar nature or late-type evolved 
stars (Paper I). While young sources are characterized by very 
steep rising continua in the infrared range, late-type evolved stars 
show more flat to falling continua. An inspection of our spectra 
shows that some of the sources seem to be composed of a falling 
SL and a rising LL spectrum. 

The two sources th at were identified as OH/IR stars by 
(ISiouwerman et al]ll99 8) both show flat spectra. The only de- 
tectable spectral features are the silicate absorption bands at 
9.8 and 18 /jm. Neither PAH emission bands nor forbidden fine 
structure lines are visible in the spectra. 

In order to identify the sources as young or late-type evolved 
objects, we looked at the slope of the spectrum as well as the de- 
tected spectral features. To allow for a quantitative classification 
of the spectral slopes, we determined the ratio of the flux densi 
ties at 33 and 5.3 fim, 



were selected in order to obtain the slope over a large part of the 
spectrum (both SL and LL parts) and to avoid spectral features at 
the short wavelengths end and the large flux jumps between data 
points at the long wavelength end. We expect large values for 
this parameter for young sources with a steep rising spectrum 
and smaller values for more evolved sources with a flat spec- 
trum. The spectra of sources that showed a flux jump between 
the SL and the LL module were first corrected for the mismatch 
and then the slope parameter was determined for these sources. 
The two OH/IR stars in our test sample have slope parameter 
values of 2 and 3. After an inspection of the data, we set the 
threshold for the identification of a source as a young object to 
be F33^/F5 3^ > 10. Sources that were only detected with the 
SLl and the LL modules or only the LL modules were identified 
by comparing their spectra to the other sources. 

In dependence of the slope parameter, we classified the 
sources as young objects or late-type evolved objects. If pos- 
sible, we subdivided the group of young objects in young stellar 
objects or HII regions in dependence of the detection of forbid- 
den fine structure lines. In the spectra of young stellar objects or 
late-type evolved stars we do not expect the detection of PAH 
emission bands or forbidden fine structure lines. In fourteen of 
the sources, the identification of the spectral features is not clear. 
However, the slope of their spectra classifies them as young ob- 
jects. 

- Young object (YO) (49 sources) 

- Young stellar object (YSO) (14 sources) 

• Steep rising spectrum 

• No PAH features or forbidden fine structure lines 

- H II region (H II) (2 1 sources) 

• Steep rising spectrum 

• Forbidden fine structure lines 

- Steep rising spectrum but no other clear features (14 
sources) 

- Late-type evolved object (LEO) (8 sources) 

- Flat spectrum 

- No PAH features or forbidden fine structure lines 

Furthermore, it should be mentioned that the classification in 
YSOs and H II regions is not completely reliable. On one hand, 
a young H II region does not have to show emission of forbidden 
fine structure lines yet. On the other hand, a young stellar object 
could have forbidden fine structure emission lines in its spec- 
trum, originating from the ambient medium or a nearby source. 
However, the slope of the spectrum unambiguously identifies 
these sources as young objects. 

Following Paper I, the flux ratio Fe/Fd of an MSX point 
source in the E (21 /im) and D (15 //m) bands can be used as an 

identification criterion for young massive stars, if Fe/^d > 2. 
For 21 of our ISOGAL sources, we could find an MSX point 
source within 8" (Only MSX sources with flux quality flags of 
2^ were considered). Based on their spectra, eighteen of them 
were identified as YSO or H II regions and their MSX counter- 
parts showed Fe/Fu flux ratios above 2. The three other sources 
were classified as late-type evolved objects, supported by the 
flux ratios of their MSX counterparts Fe/Fq of 0.85, 1.31, and 
1.57. 



The fluxes at these two wavelengths 
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Table 2: ISOGAL colors and identification of the sources. The first three columns give the source names and their [7] and [15] 
ISOGAL magnitudes. Column 4 through 6 list the [7]-[15] ISOGAL color, the spatial extent parameter cri5 and the slope parameter 
■P^. In the last two columns, the sources are identified and the detection of radio emission at their position is indicated. 



Source name 


[7] 

(mag) 


[15] 

(mag) 


[7]-[15] 


0-15 

(mag) 


F33 ljmfp53 ixm 


Identification 


Radio source? 


J 174457.4-291003 


8 


4.95 


3.05 


0.1 




Y5>0 




T1 744^0 1 OQCUi'i'^ 


f, so 


'X 40 


4 


n (\'\ 


1 "^S 








Q 

o 


A 00 


jAjo 


U.Uj 


ZjU 


vr* 
I \J 




T17A^ni A_os^^;oo 


Q 

O 


A S'i 


'X 1 7 


n 10 




II 1 1 

rl 11 




J 174504. 9-29 11 46 


4.27 


1.97 


2.3 


0.11 


58 


Hn 




J174505. 1-290937 


8 


3.83 


4.17 


0.11 


48 


YO 


— 


1174505.6-291018 


5.24 


1.78 


3.46 


0.09 


678 


Hn 




J 174506. 5-291 118 


5.29 


0.88 


4.41 


0.07 


782 


HH 




J174508.0-290655 


8 


4.34 


3.66 


0.14 




YSO 


— 


1174508.0-291039 


6.59 


4.24 


2.35 


0.11 


63 


HH 


— 


1174508.1-290840 


7.69 


5.11 


2.58 


0.06 


240 


HH 




1174508.7-290348 


8 


4.17 


3.83 


0.12 


230 


HII 


— 


J174511.3-290621 


6.44 


4.31 


2.13 


0.07 


72 


YO 


— 


J1745 12.6-290646 


8 


4.97 


3.03 


0.03 


70 


YO 


— 


J174515.4-291213 


8 


3.4 


4.6 


0.11 


1242 


YSO 


— 


J174515.9-290155 


8 


2.63 


5.37 


0.13 


246 


YSO 


— 


1174516.2-290315 


5.66 


1.03 


4.63 


0.19 


163 


YSO 


— 


1174517.1-291155 


8 


5.02 


2.98 


0.09 




YSO 


— 


1174517.8-290813 


7.3 


4.1 


3.2 


0.09 


290 


YO 


— 


J174518.1-290439 


6.3 


2.54 


3.76 


0.05 


139 


HH 




J174518.1-291051 


6.51 


3 


3.51 


0.08 


2285 


HH 


y 


J174518.3-290628 


8 


5.18 


2.82 


0.1 


22 


YSO 




J174520.0-290750 


7.34 


5.09 


2.25 


0.08 




YSO 




1174520.1-290638 


8 


4.93 


3.07 


0.12 




YO 




1174520.7-290213 


4.57 


0.4 


4.17 


0.12 


554 


YO 




1174522.1-291059 


6.81 


4.61 


2.2 


0.04 


124 


HH 




J174522.4-290242 


8 


3.68 


4.32 


0.17 


109 


HII 




J174523.3-290331 


5.52 


2.29 


3.23 


0.14 


24 


YO 




J174523.5-290225 


8 


4.62 


3.38 


0.14 




YO 




1174523.9-290310 


6.09 


3.29 


2.8 


0.14 


10 


YO 




1174524.8-290529 


8 


3.89 


4.11 


0.19 




YSO 




1174524.9-290318 


8 


3.33 


4.67 


0.14 


386 


HII 




1174525.7-290942 


8 


4.68 


3.32 


0.15 


73 


HII 




1174527.2-291126 


7.54 


5.15 


2.39 


0.1 


73 


YSO 




1174527.4-290655 


8 


5.12 


2.88 


0.13 




YSO 




J174529.5-291021 


6.89 


4.51 


2.38 


0.07 


271 


HH 




J174553.3-290406 


8 


3.55 


4.45 


0.09 




YO 




J174555.6-290521 


7.05 


3.77 


3.28 


0.09 




HH 




1174557.7-290302 


7.53 


5.05 


2.48 


0.06 


11 


YO 




1174600.1-290150 


8 


3.48 


4.52 


0.06 




HH 




1174602.5-290427 


8 


3.81 


4.19 


0.19 




YSO 




1174602.8-290359 


5.9 


2.43 


3.47 


0.08 


2 


LEO 




1174612.7-285958 


7.42 


5.08 


2.34 


0.14 


71 


HH 




1174613.1-290327 


8 


4.75 


3.25 


0.13 


2 


LEO 




1174613.5-290024 


6.84 


4.52 


2.32 


0.22 




HH 




J174617.6-285857 


8 


4.81 


3.19 


0.13 


1 


LEO 




J174618.0-290245 


8 


5.24 


2.76 


0.21 


186 


HII 





1174618.2-290136 


6.35 


0.88 


5.47 


0.19 


62 


YSO 




J174619.4-285611 


6.98 


4.57 


2.41 


0.19 


6 


LEO 




1174620.8-290532 


6.13 


4.22 


1.91 


0.04 


1 


LEO 




1174623.0-285845 


5.07 


2.97 


2.1 


0.05 


3 


LEO 




1174623.8-285300 


5.21 


3.2 


2.01 


0.03 


6 


LEO 




1174624.9-290345 


8 


5.23 


2.77 


0.09 


1 


LEO 




1174631.1-290652 


8 


5.14 


2.86 


0.09 


47 


YO 




J174631.2-285448 


6.9 


4.76 


2.14 


0.13 


17 


YO 




J174631.8-285027 


8 


4.12 


3.88 


0.21 




HH 




J174632.9-285115 


7.67 


5.21 


2.46 


0.14 


219 


HH 
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1.5 2 2.5 3 3.5 4 4.5 5 5.5 
[7] -[15] 

Figure 4: The image shows the spatial extension parameter cri5 
plotted against the [7]-[15] color The arrows indicate lower 
limits of the [7]-[15] colors for the sources which were only 
detected at 15 fim. The squares mark the sources whose [7] 
and [15] colors were not observed at the same time, diamonds 
highlight sources which were simultaneously observed at 7 and 
15 yum. The red line represents the selection criteria of Paper 
I, the blue line the selection criteria which were defined in this 
publication. The image shows that color selection criteria which 
are combined with the spatial extension parameter cri5 permit 
the distinction between young and evolved objects with a low 
contamination rate of late-type evolved objects. 



In addition, we looked for radio emission in the vicinity of 
our ISOGAL sources in images of the VLA archive at 1 .4, 4.7, 
and 8.4 GHz in order to identify potential HII regions. At the 
positions of six ISOGAL sources, which were classified as H II 
regions on the basis of their spectra, radio emission was detected. 
Several other sources are located at the edge of diffuse radio 
emission. 

Table |2] shows the ISOGAL colors as well as the spectral 
identification of the sources. Columns 1-5 list the sources names 
together with their ISOGAL magnitudes at 7 and 15 fim, the [7]- 
[15] color and the spatial extent parameter 1x15. The sixth col- 
umn contains the slope parameter pLes. jn jjjg j^st two columns 
the identification of the sources as well as the detection of radio 
emission at the positions of the sources in the VLA images are 
reported. 



2.5. Test of the ISOGAL selection criteria 

One of the intentions of this publication was to determine if the 
applied photometric color criteria ([15] < 5.25, [7]-[15] > 1.8) 
permit the selection of only young objects in a sample of 
ISOGAL point sources. The identification of eight late-type 
evolved objects in our test sample shows that the color criteria 
do not constrain the group of sources to objects of young stellar 
nature only. 

We confirmed that a combination of the spatial extension of 
the sources at 15 fim, characterized by the parameter 0-15, with 
their [7]-[15] colors gives a reliable method to distinguish be- 
tween young and late-type evolved objects. 

In Figure 21 0-15 is plotted against the [7]-[15] colors. 
Sources which were only detected at 15 fim are marked with 
an arrow. The image shows that young objects and late-type 
evolved stars cannot be divided into two separate sections in the 
plot. However, criteria which are based on the ISOGAL colors 
as well as the spatial extension parameter allow the extraction of 



young objects from a mixed sample with a small contamination 
by evolved sources. The red line represents the selection criteria 

2 < [7] - [15] < 2.5 and 0-15 > 0.1 mag 

2.5 < [7] - [15] < 3 and 0-15 > 0.05 mag 

[7] - [15] > 3, no constraint on 0-15 

which were defined in Paper I from a group of ISOGAL 
sources with known evolutionary states. Applying these criteria 
to our test sample results in the selection of 44 young objects and 
five late-type evolved stars, corresponding to a LEO contamina- 
tion rate of 10.2% (5/49). In order to reduce the LEO contami- 
nation rate by selecting more young objects and to simplify the 
selection criteria, we refined the criteria (blue line) to 

2 < [7] - [15] < 2.6 and a-15 > 0.06 mag 

2.6 < [7] - [15], no constraint on cri5 

These criteria would yield the same number of late-type 
evolved objects but 48 young objects, resulting in a LEO con- 
tamination rate of 9.4% (5/53). These criteria also agree with the 
ISOGAL source sample of Paper I and will be used in the second 
part of this publication in order to identify young sources in the 
CMZ in the ISOGAL point source catalogue. 

In order to estimate the final LEO contamination rate, also 
the remaining 11 sources for which no spectrum was assem- 
bled and the 27 sources that were identified beforehand as late 
evolved stars have to be considered. If we apply the defined 
selection criteria on the 27 rejected sources, nine of them are 
classified as young objects, increasing the contamination rate to 
22.6% (14/62). All 11 non-identified sources of our test sam- 
ple of 68 sources are sorted in the group of young objects. The 
sources labeled as group C sources are probably false detec- 
tions in the ISOGAL catalogue, increasing the contamination 
rate to 27.3% (18/66). Only one of the group D objects has an 
MSX counterpart within 8" which has a flux density ratio of 
Fe/Fd > 2. The ISOGAL sources in group D are resolved in 
more than one source at short wavelengths by the Spitzer obser- 
vations. The ISOGAL colors of these sources can be dominated 
by either young sources or more evolved objects in the cluster. 
The contamination rate is therefore 34.2% (25/73) if the group 
D objects are defined as more evolved objects and contribute to 
the contamination rate or 24.7% (18/73) if they are identified as 
young objects. In conclusion, the contamination rate with non- 
young objects is in the range 25 to 35%. 

Another spectroscopic s tudy of YSO c andidates in the CMZ 
was recently conducted by An et al.l (1201 ll) . Based on IRAC col- 
ors, they selected 107 point sources for their sample and obtained 
high-resolution Spitzer/IRS spectra for these targets. The shape 
of the absorption profile of CO2 ice at 15 //m served as their iden- 
tification criterion. In total, they identified 35 YSOs and possible 
YSOs. 

Since we only have low-resolution data for our sources, we 
could not test the An et al. identification criterion on our test 
sample. A comparison of the An et al. source sample and our 
test sample showed that both have three sources in common. 
An et al. identified only one of them as a possible YSO (source 
J174529. 5-291021) on the basis of the CO2 ice absorption pro- 
file. Beside this ice absorption band, our Spitzer/IRS spectrum of 
this source also shows emission lines of [Aril] and [Nell] in its 
spectrum. Thus, we identified this source as an H II region. We 
conclude that this source could be a transition object between 
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YSOs and HII regions, showing signatures of both phases. We 
also identified the other two sources as HII regions. An et al. 
exclude that these sources are YSOs. 

While the An et al. study focuses only on YSOs, we will 
obtain a broader picture of the star formation history in the CMZ 
since we also include objects in our test sample that are in later 
evolutionary stages of the star formation process, such as H II 
regions. 

3. Star formation in the Central IVIolecular Zone 

In this part, we will apply the defined selection criteria on all 
ISOGAL sources that have been detected in the central molecu- 
lar zone in order to determine the star formation activity in the 
central 450 pc of our Galaxy. Young object candidates are se- 
lected in an area covering +1.5° x ± 0.5° around the Galactic 
Center. The list of ISOGAL sources in the CMZ will be com- 
pared with the MSX point source catalogue (PSC) and MSX 
sources with a flux density ratio of Fe/Fq > 2 without an 
ISOGAL counterpart within 8" will be added to the CMZ source 
list of young object candidates. 

We will derive bolometric luminosities for all sources in the 
CMZ sample from their flux densities at 15 //m by using 
conversion factors which will be obtained from our test sam- 
ple sources. Furthermore, we will determine mass estimates for 
all CMZ sample sources in order to determine the total mass in 
young objects that has been formed over the last ~1 Myr and 
thus derive an average star formation rate for the CMZ. 

3.1. Selection of the Candidate Young Objects 
3.1.1. ISOGAL point sources 



Table 3: ISOGAL fields used in the present analysis. The fields 
names (Col. 1) give the galactic coordinates of their centers. 
Column 2 and 3 contain the LW filter numbers at 7 and 15 fim. 
The last two columns show the field of view (FOV) and the pixel 
size of the field. 



Name 


Filters 


FOV 




Pixel 


FC~00121-00003 


LW5 


LW9 


9' X 17' 


6" 


FC-00112-00035 


LW6 


LW9 


28' 


X 


20' 


6" 


FC-00109+00031 


LW6 


LW9 


34' 


X 


23' 


6" 


FC-00090-00003 


LW5 


LW9 


21' 


X 


17' 


6" 


FC-00062~00006 


LW5 


LW9 


12' 


X 


17' 


3" 


FC-00062-00040 


LW6 


LW9 


15' 


X 


12' 


6" 


FC-00039+00018 


LW6 


LW9 


39' 


X 


9' 


6" 


FC-00027-00006 


LW5 


LW9 


19' 


X 


17' 


3" 


FC+00004+00040 


LW6 


LW9 


46' 


X 


14' 


6" 


FC+00005-00024 


LW5 


LW9 


18' 


X 


15' 


6" 


FC+00034-00005 


LW5 


LW9 


14' 


X 


16' 


3" 


FC+00037+00017 


LW5 


LW9 


34' 


X 


10' 


6" 


FC+00059+00002 


LW5 


LW9 


14' 


X 


7' 


3" 


FC+00062-00014 


LW5 


LW9 


13' 


X 


6' 


3" 


FC-h 00066-00041 


LW6 


LW9 


46' 


X 


12' 


6" 


FC -1-00067+00038 


LW6 


LW9 


19' 


X 


14' 


6" 


FC+00089-00009 


LW6 


LW9 


18' 


X 


20' 


6" 


FC+00 124-00032 


LW6 


LW9 


23' 


X 


22' 


6" 


FC+00 127+00035 


LW6 


LW9 


23' 


X 


20' 


6" 



The data processing of the ISOGAL survey dOmont et alj 

l2003h and the creation of the ISOGAL-DENIS point source cat- 
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Figure 5: Boundaries of the ISOGAL fields used for the present 
study (see also Table O, in Galactic coordinates. Some well 
known regions are indicated: Sgr A* (asterisk), the Arches clus- 
ter (A), the Quintuplet cluster (Q), as well as the Sgr A, Sgr Bl, 
Sgr B2, Sgr C, Sgr D and Sgr E star forming regions. 



alogue have been extensively described in lSchuUer et al.l (l2003h . 
Here we will use only ISOGAL fields that have been observed at 
7 and 15 fim. A list of the fields (Col. 1) with which the sources 
were observed, the filter numbers (Col. 2, 3), field of view of the 
fields (Col. 4) and the pixel scales of the fields (Col. 5) are given 
in Table El 

The ISOGAL observations were designed to exclude strong 
IRAS sources, with F12 > 6 Jy in general, to avoid saturation 
eff'ects of the detector However, for most fields observed around 
the Galactic Center, narrow filters and 3" pixels were used, to 
allow the observations of brighter sources (up to F12 ~ 20 
Jy); nevertheless, the Galactic Center itself and a couple of other 
brightest star forming regions could not be imaged. Fig.|5]shows 
the boundaries of the ISOGAL fields, listed in Table [3] In addi- 
tion, known star forming regions are indicated. 

As shown in the first part of this publication as well as in 
Paper I, source samples based only on color selection criteria 
are contaminated with some post-main sequence stars (planetary 
nebulae, OH/IR stars with extreme mass losses) although most 
of the sources are associated with massive stars. Combining a 
color selection criterion with an indication of spatial extension 
(i.e., removing the point-like sources) allows to select infrared 
sources which are very likely to be massive (young) objects. 

The application of the selection criteria defined in the 
first part of this publication to the ISOGAL fields listed in 
Table [3] results in the selection of 485 ISOGAL sources with 
[15] < 5.25 mag. As described in subsection l2.1l and l2.3l a lower 
limit for the [7] magnitude was assigned to sources that were not 
detected at 7 fim. 



3.1.2. MSX point sources 

The Midcourse Space Experiment performed a survey of the 
Galactic plane at four mid-infrared wavelengths: 8, 12, 15 and 
21 yum, with a pixel s cale of 18" . The sensitivity of the MSX 
catalogue version 2.3 (Ega n et al]|2003h is of order 50 mJy at 
8 jum and 1-2 Jy at 15 and 21 fim. 

As has been shown in Paper I, the longest wavelengths of 
MSX provide another color criterion (Fe/Fd > 2) to select prob- 
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able young massive objects, with little contami nation by late- 
type e volved stars and planetary nebulae (see also lLumsden et aLl 
I2OO2I) . In order to have reliable flux density values for the MSX 
sources at 15 and 21 //m, we only considered sources with flux 
quality flags of 2-4 in the D and E bands. Besides the ISOGAL 
sources, we selected MSX point sources in the CMZ without an 
ISOGAL counterpart within 8" that fulfilled the specified MSX 
color criterion, resulting in 656 MSX sources. 

In total, our CMZ sample contains 1141 young object can- 
didates. Fig.|6]shows the distribution of the CMZ source sample 
in Galactic longitude I and Galactic latitude b. ISOGAL sources 
and MSX sources are marked as black and grey dots, respec- 
tively. Most of the sources in the CMZ sample are confined to 
the Galactic plane (|b| < 0.2°) and most of the MSX sources are 
located at 0° < 1 < 0.5°, which could not be observed by ISOGAL 
because of detector saturation. 

3.2. Discussion 

Among the 485 ISOGAL sources, 190 (almost 40%) have an 
MSX counterpart within a search radius of 8". The distribu- 
tions of ISO-MSX position separations (Fig. |7]) does not show 
any increase with larger separations, as would be expected from 
random associations. In addition, the agreement between ISO 
and MSX magnitudes is generally good. Thus, we do not expect 
a large contribution by false associations in our CMZ sample. 
However, the MSX magnitudes tend to be lower (the sources ap- 
pear brighter) than the ISO sources, especially at 7/8 //m. This 
can partly be explained by different central wavelengths (8 fim 
vs. 7 yum) that can make a difference for these red objects. Some 
part of the observed difference may also be due to a Malmquist 
bias: with a sensitivity more limited for MSX than ISOGAL, the 
brighter MSX sources will preferentially be associated. In addi- 
tion, the lower spatial resolution of MSX compared to ISO can 
also introduce a bias in that direction, when a single MSX source 
is resolved into several ISOGAL point sources. 

Among these 190 MSX sources, 177 have flux quality flags 
of > 1 in A and C band and > 2 in D and E band. Their positions 
in the [15]-[21] vs. [8]-[15] color-color diagram are shown in 
Fig. [H The distribution of all MSX sources in the region de- 
tected in the 4 MSX bands (1833 sources) is shown as small 
black dots whereas the 167 sources with an ISOGAL counter- 
part within 8" are shown as black boxes. The 656 MSX sources 
with Fe/Fd > 2 without an ISOGAL counterpart within 8" are 
marked with grey circles. As seen in Fig. [8] almost half of the 
MSX sources in the CMZ appear in a 'clump', around [8]- 
[15] ~ 2 and [15]-[21] ~ 0.5. Only ~ 20 of our ISO-MSX 
sources are located in this clump, and could actually be 'nor- 
mal' field objects. On the opposite, 80% of the ISO-MSX ob- 
jects have a Fe/Fd ratio above 2. Therefore, we conclude that our 
candidate young massive objects, selected on the basis of their 
ISO magnitudes, show peculiar MSX colors, and indeed trace a 
specific population. 

3.2.1. Interstellar extinction 

Interstellar extinction in the direction of the CMZ can exceed 
30 mag in the V band. Since we want to estimate the bolometric 
luminosities of our sources from their flux densities at 15 yum, 
we have to correct them for interstellar reddening for all sources 
in the CMZ sample. 

We used an extinction map, computed by Schulthe is et alj 
(I2009h . which is based on the combination of photometric data 
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Figure 8: MSX Color-Color diagram. All MSX sources with flux 
quality flags > 1 in A and C band and > 2 in D and E band 
in the CMZ are plotted with black dots. The black boxes mark 
afl MSX sources with an ISOGAL counterpart within 8". MSX 
sources with Fe/Fu > 2 without an ISOGAL counterpart within 
8" are shown as grey circles. The image shows that almost half 
of the MSX sources in the CMZ appear in a 'clump', around 
[8]-[15] ~ 2 and [15]-[21] ~ 0.5, but only ~ 20 of the MSX 
sources with an ISOGAL counterpart are located in this clump, 
and could actually be 'normal' field objects. On the opposite, 
80% of the MSX objects with an ISOGAL counterpart have a 
Fe/Fu ratio above 2 and indeed trace a specific population. 



from the 2MASS and IRAC/GLIMPSE surveys. We derived a 
value for Ay toward each object in our CMZ source list as an 
average of all data points within a radius of 0.5' of the ob- 
ject position. For CMZ sources that are not located within the 
area of this extinction map, we used an older extinction map 
dSchultheis et al.l[T999h which is based on photometric data from 
the DENIS survey. 

Furthermore, we extracted -r^ conversion factors from the 

results of lJiang et al.l (l2005l) who give a value for each ISOGAL 
field. For those MSX sources that are not located in an ISOGAL 
field, an average value of 0.4 was set as the^ conversion fac- 

tor. The -r^ conversion factor was obtained from the results of 

iNishivama et~an (l2008h (^ = 0.062). This value was derived 
specifically for the Galactic Center environment. 

The corrected magnitudes at 15 yum of the ISOGAL sources 
were then converted into flux densities. The MSX flux density 
values were corrected by translating the flux density into magni- 
tude values, subtracting A 15 and then retranslating the corrected 
magnitude values into flux densities. 

3.2.2. Bolometric luminosities 

In order to derive a conversion factor between the flux density 
Fi5 and the luminosity Lboh we fitted the Spitzer/IRS spectra of 
the 49 y oung objects in our test sample with the SED fitting tool 
of Robit aille etaP (12007 ). This tool contains 20000 YSO mod- 
els, covering a wavelength range of 0.01 to 5000 jum, calculated 
at 10 different viewing angles and 50 circular apertures, which 
results in a total number of 10 million SEDs (Robitaille et all 
2006). The models cover a large stellar mass range of 0.1 to 50 
M0 and all evolutionary stages from pre-stellar to pre-main se- 
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Figure 6: Areal coverage of our CMZ sample with Galactic longitude 1 and Galactic latitude b. Black dots mark the ISOGAL 
sources, grey dots the MSX sources without ISOGAL counterparts. The image shows that most of the MSX sources are located at 
0° < 1 < 0.5°, which could not be observed by ISOGAL because of detector saturation. Most of the sources in the CMZ sample are 
confined to the Galactic plane (|b| < 0.2°). 
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Figure 7: From top-left to bottom-right, distributions of: Angular separation between the ISOGAL and the associated MSX source, 
ISO-MSX 7 jjm magnitude differences, ISO-MSX 15 /im magnitude diff'erences, and ISO-MSX differences in colors. 



quence. During the fitting, the tool applies all the models to the 
entered data, leaving the distance and the external foreground ex- 
tinction Ay as free parameters within a user specified range. The 
result is a list of 10000 fits that are characterized by ax^ value. 
Since the parameter space is sampled too sparsely to allow an 
accurate determination of the surface minima and, therefore. 



the appropriate confidence intervals, a specified criterion for 
"good" fits has to be selec ted. We defined the fo llowing crite- 
rion: -xi:,i < 3 ■ «points (iRobitaille et a l.l l2007h withxl,,i be- 
ing the x^ value of the best fitting model and npoints the number 
of data points used in the fitting process. Although this selection 
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criterion is kind of arbitrary, a stricter criterion would imply an 
over-interpretation of the fitting results. 

Since the distance and the external foreground extinction 
serve as free parameters, a distance range and a range for the 
visual extinction Ay have to be specified. We assumed that all 
CM Z sources are loc ated at the Galactic Center distance (8.0 
kpc. lReid et al.ll2009l) . Accordingly, we chose the fitted distance 
range to be 7-9 kpc. We set Ay to be 15 to 50 mag, corre- 
sponding to typical maximal and minimal values of Ay in the 
area of the 49 sources, extracted from the extinction map of 
ISchultheis etal] (|2009). 

For the fitting, we selected between six and eleven data 
points per spectrum, spread between 5 and 35 fim, to well reflect 
the underlying continuum. Since the SED tool does not model 
emission lines, PAH or ice features, we avoided the selection of 
data points in thus affected regions of the spectrum. Since the 
fitting tool permits to enter data taken with different apertures 
we can select the data points from the spectrum without correct- 
ing for the flux jumps that are visible in some of the spectra. For 
the size of the apertures, we chose the widths of the SL and LL 
modules. 

We selected the models following the above defined criterion 
and calculated the bolometric flux densities using the luminosity 
values and the distance values of the different models. Then, we 
reconverted them into luminosity values, using a single distance 
of 8.0 kpc for all models. We calculated the weighted average 
and standard deviation of the bolometric luminosities for each 
source with the inverse of the;^^ values as the weights. 

In the next step, we plotted the obtained bolometric lu- 
minosities of the 49 ISOGAL sources against their cor- 
rected ISOGAL flux densities at 15 /urn (left image in 
Fig. |9]l. The standard deviations served as error bar val- 
ues. A linear fit to the data provided the conversion factor 
(§;7)isoGAL = (4447 + 106) Lo/Jy. For the 18 ISOGAL sources 
with an MSX counterpart within 8" we plotted their bolomet- 
ric luminosities against the corrected MSX flux densities at 
15 fj.m (right image in Fig. |9l). The slope of the fitted line is 
(^)msx = (3736 + 209) Lo/Jy. 

Afterwards, we calculated the bolometric luminosities for all 
sources in our CMZ sample from their Fis flux densities. The 
luminosity values range from 0.5 ■ lO-' to ~ 10^ Lq, with a total 
of -2.5 ■ 10^ Lq. 

3.2.3. Masses 

It is commonly accepted that hydrogen burning in high-mass 
stars sets in while the star is still accreting material from the 
circumstellar birth cloud. At this phase, the surrounding mate- 
rial is heated and eventually ionized by the UV radiation of the 
star. Based on our ISOGAL selection criteria, the sources in our 
CMZ sample are selected to have bright mid-infrared emission 
from warm dust and hence are still highly embedded in their 
birth clouds but already contain a central heating source. We can 
therefore assume that most of the sources in our CMZ sample are 
likely to be intermediate- to high-mass stars which have already 
started hydrogen burning and thus are akeady on the zero age 
main sequence (ZAMS). Assuming that the bolometric luminos- 
ity of each source is that of a single ZAMS star, mass estimates 
for all CMZ sample sources have been derived, using the relation 
(e.g. Eddington. 1926; .Grifliths et aLil988,: .Bhattacharva et alj 
12001 

Lq ~^Mq' ■ 
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Figure 10: The image shows the distribution of the source masses 
in our CMZ sample. The line represents the Kroupa IMF 
^(M) - ^oj(M) ■ M"^'^. The histogram clearly shows that low- 
mass and very massive objects are missing in our CMZ sample 
due to our selection criterion [15] < 5.25 mag and saturation ef- 
fects, respectively. 

The distribution of the masses, which ranges from 6 to 
58 Mq, is roughly consistent with a Salpeter or Kroupa initial 
mass function (IMF.IS alDeteJT95llKroupal200TI) . However, the 
histogram (see Fig. [TOl i clearly shows that low mass objects are 
missing from our CMZ sample which is due to the selection cri- 
terion [15] < 5.25 mag, corresponding to a lower mass limit of 
6.3 Mq. While most of the total luminosity in a young popula- 
tion comes from high mass stars, most of the total mass is con- 
tained in the high number of low mass stars. Since we are only 
selecting the most massive objects (the most luminous ones), we 
miss the stars with lower masses. In order to obtain a complete 
picture of the mass distribution in th e CMZ, we wi ll extrapolate 
the IMF to lower masses, following (lKroupal2001h 

^(M) = ^0,1 (M)M"--^ for O.5M0 <M < IIOMq 

and 

^(M) = ^o,2(M)M"'-^ for O.O8M0 <M < 0.5Mq. 
and 

^(M) = ^o,3(M)M"°-^ for O.O1M0 <M < O.O8M0. 

We fitted the histogram with a curve of the form 
y = ^o,i(M) ■ x-2-^ and obtained ^o,i(M) = 16019 M^^ 
(black line in Fig. [TOl l. In order to obtain a continous IMF, 
we derived ^o,2(A^) from ^q,\{M) at M = 0.5 M0 and ob- 
tained ^o,2(M)' = 32038 MP^ and ^(),3(M) from ^o,2(M) at 

M = 0.08 M© and obtained ^o^M) = 400475 M^^'-^ The to- 
tal mass of young objects in the CMZ can then be derived with 
Mtot = / M ^(M) dM. In the mass range 0.01-120 Mg, the 
total mass is ~ 77000 M.q. 

The histogram can roughly be divided into three mass 
ranges: 0-8 Mq, 8-20 Mq and 20 - 60 Mq. As explained 
above, low-mass objects are missing in the lower mass range 
due to our selection criterion [15] < 5.25 mag. In the second 
mass range, the number of sources is underestimated by the fitted 
IMF by roughly 35% compared to the number of sources in our 
CMZ sample. Confusion due to the high population in the CMZ, 
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Figure 9: The images show the relation between the bolometric luminosity and the extinction corrected ISOGAL (left) and MSX 
(right) flux densities at 15 yum for the young objects in our test sample. The error bars correspond to the standard deviations of 
the bolometric luminosity values. The slopes of the fitted lines are (4447 + 106) L0/Jy (left) and (3736 + 209) L0/Jy (right), 
respectively. 



saturation effects of the detectors as well as rareness limit the de- 
tection of very luminous, and therefore massive, sources, which 
would explain the lack of objects in the upper mass range. As ex- 
plained in Part 1 of this publication, we expect a contamination of 
our ISOGAL-CMZ sample with late-type evolved objects with a 
rate of 25% - 35%. If we assume that the MSX-CMZ sample is 
also contaminated with more evolved objects with a similar rate, 
we estimate the total contamination rate of our CMZ sample to 
be roughly 30%. Thus, we expect to have ~ 30% more sources 
in our CMZ sample over the whole mass range than predicted by 
the fitted IMF. This effect is not noticeable in the first and third 
mass range due to the number of missing objects. However, in 
the second mass range the larger number of sources in our sam- 
ple compared to the number of sources predicted by the fitted 
IMF can be explained with the contamination of our sample with 
late-type evolved objects if we assume that our sample is com- 
plete in this mass range. 



3.2.4. Stellar clusters and the Sagittarius star forming 
regions In the CMZ 

Several regions in the CMZ are known as active star forming 
sites such as the Sagittarius star forming regions (Sgr B2+B1, 
Sgr C, Sgr D, Sgr E) and the three massive stellar clusters that are 
close to the Galactic Center (Arches, Quintuplet, Central clus- 
ter). To obtain an estimate of the number of sources in our CMZ 
sample that are located in these active star forming regions and 
therefore the influence of these sites on the star formation history 
in the CMZ, we counted the number of ISOGAL sources and 
MSX sources without an ISOGAL counterpart within the area of 
each of these star forming regions. In Sgr B2 and Bl we found 
35 ISOGAL and 32 MSX sources with a total mass of 930 Mq. 
We counted 6 ISOGAL and 6 MSX sources in Sgr C with a total 
mass of 201 M0. In Sgr D, we detected only one MSX source 
with a mass of 42 M0. Neither ISOGAL nor MSX sources were 
detected in Sgr E as well as in the three stellar clusters. These re- 
gions are probably more evolved than our CMZ sample sources 
since they are several million years old and their young stars are 
not anymore deeply embedded in their natal cloud which makes 



them appear less red. Thus, our sample traces more recent star 
formation activity. 

3.2.5. Average star formation rate 

In order to appear bright at mid-infrared wavelengths, the 
sources in our CMZ sample still have to be deeply embedded 
in the dust of their birth clouds. Following IWood & Churchwelll 
(Il989l) . this phase lasts about 10% of the lifetime of an O-type 
or early B-type star, which corresponds roughly to ~ 1 Myr for 
a typical BO star Assuming that all objects in our CMZ sample 
have formed over an average time span of 1 Myr, we can convert 
the total mass of young objects (~ 77000 Mo) to an average star 

formation rate of 0.08 M0 yr_^ 

Previous studies ( Croc ker et al. ] I20T0I: lYusef-Zadeh etaP 
2009) published a range of values fo r the star forma tion rate 
in the CMZ from 0.08-0.15 M0 vr '. lAn et all dloill) make a 
rough estimate of the star formation rate based on their sources 
which are in common with the Yusef-Zadeh et al. source sam- 
ple, deriving a value of 0.07 M0 yr"' . Although our calculations 
are based on a number of simplifying assumptions (e.g. linear 
relation between F15 and Lboi, transformation of the bolometric 
luminosities into masses of single ZAMS stars, ...), we can state 
that our result is consistent with other independent studies of the 
star formation rate. 

4. Summary 

The goal of this study was to test and revise ISOGAL selection 
criteria on a test sample of 68 unknown ISOGAL sources and to 
apply these criteria to all ISOGAL point sources in the CMZ in 
order to find young objects and derive the average star formation 
rate over the last 1 Myr We obtained Spitzer/IRS observations 
for a test sample of 68 ISOGAL sources and assembled spectra 
in the wavelength range 5-38 fim for 57 of these sources. In 
dependence of the detection of PAH or forbidden fine structure 
line emission and the slope of the spectrum, we classified the 
sources as young or late-type evolved objects. 

After revising the ISOGAL selection criteria, which are 
based on the ISOGAL [7]-[15] color and the spatial extent pa- 
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rameter cri5, we applied them to all ISOGAL sources in the CMZ 
(± 1.5° X ± 0.5° around the Galactic Center), selecting 485 
sources. Furthermore, we added 656 MSX sources in the CMZ 
which fulfilled the criterion Fe/Fd > 2 (Fd, Fe flux densities in 
band D (15 fim) and E (21 fim)). 

After fitting the SEDs of the 47 young objects in our test 
sample, we obtained bolometric luminosity values and plot- 
ted them separately against the extinction-corrected ISOGAL 
and MSX flux densities at 15 yum of the sources. We obtained 
^ conversion factors and applied them to our CMZ sample. 
Assuming that the bolometric luminosity of each source cor- 
responds to the luminosity of a ZAMS star, we calculated the 
mass for each source. The distribution of the masses roughly 
follows a Salpeter or Kroupa IMF but low-mass and very mas- 
sive objects are missing from the CMZ sample. The total mass of 
young objects in the CMZ (~ 77000 Mq) was calculated from 
a Kroupa IMF between 0.01 and 120 M0. Since the detection 
as bright mid-infrared sources requires that the sources are stiU 
deeply embedded in their dust cocoon, the sources still have to 
be very young. Assuming that this total mass in young objects 
has formed over the last 1 Myr, we obtained a star formation rate 
of 0.08 M0 yr"^ for the CMZ which is consistent with previous 
studies. 
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Figure 3: Background subtracted spectra. The data of the different subslits are presented by the three different colors (green: SL2 or 
LL2, blue: SL3 or LL3, red: SLl or LLl). 
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Figure 3: Continued 



J1 74508.0-291 039 



J1 74508.1 -290840 



1000 
100 
10 
1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
Order 1 




10 15 20 25 30 35 40 

X [|im] 



1000 
100 
10 

1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
" Order 1 













10 15 20 25 30 35 40 

X [urn] 



J1 74508.7-290348 



J1 7451 1.3-290621 



1000 
100 

10 
1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
Order 1 




10 15 20 25 30 35 40 

X [|im] 



1000 
100 
10 

1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
Order 1 




10 15 20 25 30 35 40 

X [urn] 



J1 7451 2.6-290646 



J1 7451 5.4-291 21 3 



1000 
100 

10 
1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
Order 1 



10 15 20 25 30 35 40 

\ [|im] 



1000 
100 
10 

1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
Order 1 




10 15 20 25 30 35 40 

A, [urn] 



J1 7451 5.9-2901 55 



J1 7451 6.2-29031 5 



1000 
100 

10 
1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
Order 1 





10 15 20 25 30 35 40 

X [|im] 



1000 
100 
10 

1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
Order 1 



10 15 20 25 30 35 40 

X [urn] 



K. Immer et al.: Recent star formation in the inner Galactic Bulge II - The CMZ, Online Material p 3 

Figure 3: Continued 
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Figure 3: Continued 



J1 74620.8-290532 



J1 74623.0-285845 



1000 
100 
10 
1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
Order 1 



10 15 20 25 30 35 40 

X [|im] 



1000 
100 
10 

1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
Order 1 




10 15 20 25 30 35 40 

X [urn] 



J1 74623.8-285300 



J1 74624.9-290345 



1000 
100 

10 
1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
Order 1 




w 


/ \^ 



10 15 20 25 30 35 40 

X [|im] 



1000 
100 
10 

1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
Order 1 



10 15 20 25 30 35 40 

X [urn] 



J1 74631. 1-290652 



J1 74631. 2-285448 



1000 
100 

10 
1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
Order 1 




• ^ 

- ^ ■ '. .;' 





10 15 20 25 30 35 40 

X [|im] 



1000 
100 
10 

1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
Order 1 



10 15 20 25 30 35 40 

A, [urn] 



J1 74631 .8-285027 



J1 74632.9-2851 15 



1000 
100 

10 
1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
Order 1 



1000 
100 
10 

1 

0.1 
0.01 
0.001 



Order 2 
Order 3 
" Order 1 




_ ' * * "**, • V 


k 



10 15 20 25 30 35 40 

X [|im] 



10 15 20 25 30 35 40 

X [urn] 



